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Abstract 

Canonical quantization of electromagnetic field inside the time-spatially dis- 
persive inhomogeneous dielectrics is presented. Interacting electromagnetic 
and matter excitation fields create the closed system, Hamiltonian of which 
may be diagonalized by generalized polariton transformation. Resulting dis- 
persion relations coincide with the classical ones obtained by the solution of 
wave equation, the corresponding mode decomposition is, however, orthogo- 
nal and complete in the enlarged Hilbert space. 
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I. INTRODUCTION 



The investigation of electromagnetic field in dielectrics has attracted growing attention 
recently [|l]-[5|]. Particularly, quantum aspects of this problem are of current interest due to 
the potential applications in technology of nanostructures. This research includes investi- 
gation of quantum wells embedded in microcavities |J and generation and propagation of 
nonclassical states of light. In this connection, the quantization of electromagnetic field in 
inhomogenous time-spatially dispersive linear medium represents a nontrivial problem 0. 

Our considerations are motivated by standard electromagnetic theory ||. Suppose for 
concreteness the geometry of closed cavity (R=l) with dispersive inhomogeneity (refractive 
index n(fl)) along the z-axis as sketched in the Fig. [I]. For the sake of simplicity the 
s-polarization of electric field only will be assumed in the folloving. Using the Maxwell 
equations, the eigenmodes of the cavity with the time dependance e tQt may be specified as 
solution of Helmholtz (time independent) wave equation 
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This equation may be simply solved in the regions where the coefficients are continuous 
function of z. The electric field E must be continuous together with its first derivation |jj 
for each \z\ < L/2. The boundary conditions at ±L/2 are given as E(±L/2) = and they 
yield the dispersion relation 
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where 



Ql = %~q 2 and Q' 2 = n 2 (9)% - q 2 . 



Here q represents the 2D component of wave vector paralel to the boundaries. This tran- 
scendent equation may be solved yielding the discrete set of eigenvalues. Nevertheless, a 
simple analysis shows that the corresponding eigenfunctions E q ^ m (z) are not orthogonal since 
the "potential" x(^) depends on the frequency. This is the source of theoretical troubles, 
since the decomposition of electric field is becoming questionable. The quasinormal modes 
in leaky macro cavity were used in Refs. P~P|- 

The purpose of this contribution is to clarify the quantum meaning of this classical elec- 
tromagnetic problem. Adopting the canonical quantization scheme formulated by Huttner 
and Barnett ||, we will formulate and solve the problem of quantization of electromag- 
netic field in closed cavity with dispersive inhomogeneity. Normal modes are associated 
with (generalized) polariton transformation. Two extreme cases of this formulation may 
be distinguished as problem of exciton confined in quantum well embedded in microcav- 
ity || and the above mentioned classical problem of leaky macrocavity. Polariton solution 
exactly yields the orthogonal decomposition. Additional degrees of freedom are associated 
with matter excitations. For the sake of simplicity all the quantum considerations will be 
performed for special form of singular refractive index. In Appendix |A] this result will be 
extended to the general form of an arbitrary refractive index fulfilling the Kramers-Kronig 
relations. 

II. CANONICAL QUANTIZATION 

Let us formulate canonical description of interaction of transversal electromagnetic field 
with matter. Neglecting other losses the Lagrangian reads 

L = Jd 3 r £(r), (3) 

where Lagrangian density is £ = C ern + C mat + C int , 

£ em = |[A 2 -c 2 (Vx A) 2 ], (4) 

£-mat — 7^ [X 2 — CUqX 2 ], (5) 
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boldface characters denote vectors and dot means time derivation J^. Electromagnetic part 
is represented by vector field A defined in the whole cavity. Polarisation part is modeled 
by harmonic oscillator field with amplitude vector X, which is non-zero only in the interval 
of inhomogeneity \z\ < 1/2, ujq being the frequency and p density. Interaction of both fields 
is characterized by interaction constant a. For the simplicity, linearly polarized fields with 
polarization paralel to discontinuity planes z = ±1/2 will be assumed (s-polarization) in the 
following. The vector field A may be interpreted as electric intensity and both fields may 
be represented as (real) scalar fields. The Lagrange- Euler equations then read 
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yielding exactly the wave equation ([]]) for susceptibility 
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A. Free field decompositions 

The interacting fields may be quantized using expansion in orthogonal basis relevant to 
respective free fields. Standard approach may be used in the plane xy perpendicular to 
the direction of inhomogeneity, since q is conserved due to the translation symmetry. In the 
^-direction, the eigenfunction for light and matter excitation parts should be distinguished. 
Assuming e* qT /(27r) dependence, r being the projection of 3D r vector into the x?/-plane, 
{ip m (z)} are the solutions of time-independent wave equation 
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? m (z) = 0; Ql = ^f- q \ (10) 



fulfilling the given boundary conditions on the interval \z\ < L/2. Assuming for concreteness 
perfect reflection on the end mirrors, we have Q m = irm/L; m = 1,2,3, .... Consequently 



frequencies are quantized as £l m (q) = cJ (rem/ L) 2 + q 2 , q = |q|. Corresponding eigenfunc- 
tions are given as ip m (z) = a/2 sm[mii(z/L + 1/2)]. The inhomogeneity of matter excitations 
is included in the definition of eigenf unctions x q ^(r), since they are non-zero only on the 
interval \z\ < 1/2. Two sets of functions <^ q , m (r) and Xq,?( r ) defined in the 3D space are 
orthogonal and complete in the volumes of quantization 

/ d 3 r ^ q>m (r)^*, n (r) = 6(q - q')5m,n, E / rf2( l < ^q > m(r)<^ q , m (r / ) = 5(r - r'), 
J dh Xq, 5 (r)x q ,,(r) = *(q - q')<%„ £ / d\ X^)x^') = S(t - v')B{z). 
The cross-products are given by matrix elements 

/ V^q.ml^Xq'^ (r)c? 3 r = <J(q - q) i^ m , f , 



(11) 

Here the full set of functions {x^i z )} create an orthogonal and complete system of functions 
on the interval \z\ < 1/2. The explicit form of correct exciton wave functions Xq,t( z ) re P" 
resents very complex problem, solution of which is beyond the scope of this contribution. 
The basis used in decomposition of matter excitations will be considered as given and index 
£ will be intuitively interpreted as energy levels of free matter excitations. In the following 
the wave functions of electromagnetic field will be consistently enumerated by Latin indices, 
whereas Greek ones will be used for the matter excitations. 

Even if we started with the Lagrangian of harmonic oscillator field (|^), we will incorporate 
into our scheme also a little more general models corresponding to the finite number of 
energy levels £ in the decomposition of harmonic oscillator. This technique corresponds 
to quantization with confinements tending to spatial dispersion. Particularly, the matter 
excitations in the ground state only correspond to the case of single exciton confined in 
quantum well 0. For this purpose, the range of sumation in respective decompositions of 
matter excitations will not be specified explicitly and will be mentioned in the discussion of 
the final results only. The respective expansions used in the following then read 
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B. Hamiltonian formalism 



The total Lagrangian L = L em + L mat + L int may be quantized as 
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Here the prime means the integration over the half of the reciprocal space. Canonically 
conjugated variables are given as 
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and as the complex conjugated relations. Hamiltonian reads 
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Standard quantization is prescribed by commutation relations between operators 
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^q^qV; =^%^5(q-q'). 
Annihilation operators of electromagnetic field 
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and matter excitations 



2%UJn 



are fulfilling the ordinary boson commutation relations 
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The definition may be extended to full space of q vectors as 
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Hamiltonian then reads 



H = H / rf2 q hVl m (q) a qjm a q , m + ^ / d 2 q ftu; & qif 6 q ,g + 



m 

%h 



K r , 



+ ^-GV^oX! /" cPq ^i- (a q , w + al qjn ) • (6 q>5 - 6_ q , 5 ) 



+ ^ 2 £ / A 



fi m (g)fj n (g) 



v a q,n + a -q,n) ' ( a q,m + a — q,m)) 



where 



(25) 



(26) 



(27) 



(28) 



and the effective interaction constant is abbreviated as 



(29) 



G = a/y/e^p. 



(30) 



The Hamiltonian of this type has already been investigated recently in connection with 
polariton effects || and exactly corresponds to the many-exciton lines interacting with 
eletromagnetic field discussed in Ref. E 
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III. DISPERSION RELATIONS 



General form of polariton transformation diagonalizing the Hamiltonian is given as 



5 q , n = Yl + J2 X ^K,^ + Y^ma-^m + Yl Z ^ b -<t,f ( 31 ) 

m £ m £ 

Index m exhausts all the cavity modes and £ similarly does all the modes of decomposition 
of matter excitations. Standard diagonalization condition 



TiVLB, 



(32) 



yields the dispersion relation for eigenfrequency Q and relations for coefficients in ([H]). The 
anticipated operator solution is normalized with respect to the boson commutation relation 



(33) 



Straightforward but lengthy calculations lead to the following equations (dependence on q 
will be omitted for brevity) 

(Q m - Q)W m + l -G 2 £ ~^={W k - Y k ) + ^G^L £ K m ^ + Z 6 ) = 0, (34) 
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(n m + n)Y m - ^g 2 £ wPfr ( w * ~ Y ^ ~ ^ G ^r £ + z d = o, (35) 
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(Q - cu )X ( + % -G^Y. - Y m ) = 0, (36) 



(n + uo)Zt - ^V^]T %^(W m - Y m ) = 0. (37) 

Finally, the dispersion relation follows as the condition for existence of non-trivial solution 
of linear equations for T m 

G 2 tt 2 
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where T m = (W m - Y m )/y/U^ 
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The algebraic treatment may be equivalently replaced by an analytical method. The 
recurent system of linear equations fl5S|) may be rewritten to the form of integro-differential 
equation. Hence the quantum problem is related to the classical solution of wave equation 
with dispersive inhomogeneous medium. Let us define formally the function 

which is continuous and has continuous derivation dA/dz on the interval \z\ < L/2. After 
simple manipulations using relations fllOD , ( |Tl~D and ( |29|) we find that it fulfills the equation 



d 2 2 1 r<-<^ 



A (z) = -»(-) 2 , g m ! ' ExsMxsM^M dz' , (39) 
c yujQ — si j J— 1/2 £ 

where Q 2 Z = (fi/c) 2 — q 2 . This alternative representation of (|38|) may be interpreted as scalar 
wave equation with the time-spatially dispersive inhomogeneity along the z-axis. Conse- 
quently, we are able to associate the quantum problem-diagonalization of the Hamiltonian 
(f28|), with the classical solution in electromagnetic theory. This analogy represents a power- 
ful tool for specification of dispersion relations (TT3| . The wave equation may be solved 
separately in the regions where all the functions are continuous. Dispersion relation is then 
given as necessary condition for continuity of the solution and its first derivation on the 
boundary \z\ = 1/2. The general solution of eq. (|39"D on the interval \z\ < 1/2 is given as 
superposition of particular and fundamental solutions 

A(z) = E c ? / dz ' G ( z ' z ')xd z ') + A 2 e iQ * z + B 2 e~ iQ * z , (40) 

A2 and B2 being general multiplicators in fundamental solution, G(z, z') being the Green 
function of the operator d 2 /dz 2 + Q 2 . We may take the explicit form 

G(z, z') = -— sm(Q z \z - z'\). 

Coefficients of particular solution are given in accordance with the relations 

G 2 Q 2 f G 2 Vt 2 ( f 

^ = ~ c 2 (u 2 -n 2 ) J dz *d z ) A ( z ) = - c2(cu 2_^ 2) {E c ^ J dz dz'xt(z)G(z,z') X ^z') 



(41) 

+A 2 J dz' xd z 'y QzZ ' + B 2j dz' xd z ') e ~ iQzZ '}- 



Parametrizing the solutions on the intervals L/2 < \z\ < 1/2 as Aje l ® z + Bje~ l ® z ; j = 1,3 
the dispersion relation may be found in closed form for limited number of terms £ for any 
mode decomposition x^( z )- 

Algebraic method will be used in the following sections [111 A| and [111 B , whereas the 



analytical method will be applied on the analogy of classical problem in the section |111 lj . 

A. One exciton dispersion relation 

Closed form of dispersion relation may be find easily in some special cases. Particularly, 
exciton confined in quantum well (QW) is characterized by the only energy level in the 
expansion of matter excitations £ = (ground state of excitations). The integro-differential 
equation ( |3~9"1) indicates the spatial dispersion. Nevertheless, the dispersion relation may 
be found without solving it, since the coefficients D m ^ are factorized as product D njTn = 



K n fl ■ K m: o. Equations (|38|) then directly yield the necessary condition 



G 2 Vt 2 ^-^ Kmfi 



^o 2 -fi 2 ?^-fi 2 (42) 
representing the dispersion relation of QW polaritons embedded in microcavity 0. 

B. Many exciton dispersion relations 

Dispersion relation for many excitons interacting with electromagnetic field will be 
demonstrated on the explicit example of two excitons (£ = 0,1). Taking into account 
the form of the kernel D m n = K mt0 K nt0 + K m iK n l , the system of linear equations ( EE ) may 
be rewritten as 

^ G 2 Vt 2 , -, Km,Q . G 2 Vl 2 x K m fiK m> i 

x ° ~ ^l^w e m -n 2 + Xl ^a 2 ~ e m-n 2 



_ G 2 Vl 2 x K m flK my \ G 2 Vt 2 , ^m,i 

where x% = ^Z m T m K m> ^ £ = 0, 1. Dispersion relation for two exciton then reads 
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The sums involved in the dispersion relation may further be expressed in the analytical form 
as done in Ref. |J, but this is beyond the scope of this paper. Let us only note that the 
application of analytical method yields directly the closed form of dispersion relations. 



C. Classical dispersion relations 

As the last special example, the dispersion relations resulting from the classical electro- 
dynamics in dispersive inhomogeneity (|I|) will be considered from the quantum viewpoint. 
This case is characterized by the decomposition of matter excitations, which is complete on 
the interval of inhomogenity (i.e. all the energy lines of excitations are included). Relation 
( ^9f) then reads 

rl/2 rl/2 rl/2 

Ai,m = E / ( Pn(z)xd Z ) dz ' / ( Pm(z f )Xt( z ') dz ' = / <Pn(z)<Pm(z)dz. (45) 
^ J -1/2 J -1/2 J —1/2 

Spatial dispersion in ( j39|) disappears yielding scalar wave equation identical with Lagrange- 
Euler equation (0) for the susceptibility @ and for the dispersion relation (||). Consequently, 
the classical solution yielding non-orthogonal eigenmode functions was completed by the 
fully quantum treatment characterized by the same dispersion relations but orthogonal de- 
composition. 



IV. CONCLUSION 

The problem of canonical quantization of electromagnetic field in linear dispersive in- 
homogeneous media was formulated in terms of overlaping of wave functions related to 
quantization of free electromagnetic and matter excitation fields. Diagonalization is given 
by generalized polariton (Hopfield) transformation. Within the classical electrodynamics, 
it may be also interpreted as the wave in spatially dispersive medium. The description of 
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macroscopic inhomogeneity and quantum well polaritons may be unified in this way into the 
same framework. 

Even if quantum and classical problems yield the same dispersion relations, there is a dif- 
ference between both the treatments. Since the electromagnetic field itself is not conserved, 
the respective eigenfunctions are not orthogonal representing the system of quasinormal 
modes P-pl- They may be used for description of electromagnetic field inside the cavity, 
however since completeness and orthogonality relations should be redefined, the description 
is more complicated than in the ordinary case of orthogonal modes. On the other hand the 
quantum solution is characterized by an ordinary orthogonal diagonalization in the form 
of generalized polariton transformation acting on the space of coupled electromagnetic and 
matter excitation fields. Standard description of the time evolution may be used, since the 
normal modes are orthogonal and complete. 
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APPENDIX A: INHOMOGENEITY WITH AN ARBITRARY REFRACTIVE 

INDEX 

The theory developed above describes the lossless dispersive inhomogeneity characterized 
by the real part of the susceptibility (|) 

= ^ < A1 > 

Principle of superposition may be used now to get an arbitrary susceptibility as done in 
Ref. |14j| . The electromagnetic field will interact with an ansamble of independent oscilla- 
tors with Lagrangian densities as in (HH), distinguished by frequencies u ,Ui,.... etc. and 
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by different parameters ao,«i, etc. and po, Pi, etc.. In the limit of continuous frequency 
distribution u, the parameters are considered as frequency dependent a(to), p(to) yielding 
the total contribution to the Lagrangian density 

C mat = Cd^in-^ll (A2) 



o 



oo 



tint = - dua{uj)A ■ X w . (A3) 

JO 

Further development of the theory runs similarly as in the single-frequency case: Matter 
operators are denoted by modal index £ and by an additional continuous index uj. The 
effective interaction parameter fl30[) is freqency dependent G(u). The final system of linear 
equations (p8|) reads 

O 2 roc O^d A 

what corresponds to the real part of the susceptibility 

r°° G 2 (uj) 

*™=L dw ra- < A5 > 

Due to the Kramers-Kronig relations jTJ]], the real and imaginary part of susceptibility 
are mutually related by Hilbert transformation as 

x ^) = ir^U, ( A6) 



IT JO LO 



IT JO LO 



The imaginary part of lossless and general cases therefore read 

X "(n) = -ixG 2 5{VL 2 - ul) (A8) 

and 

X"(fi) = -7T / duG 2 (co)S{uj 2 - Q 2 ), (A9) 
Jo 

respectively. The relations ( |A5| ) and ( |A9| ) represent the desired extension of the theory 
with single-frequency of matter excitations and singular susceptibility into an arbitrary 
case characterized by general refractive index (susceptibility). 
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FIGURES 




FIG. 1. Geometry of closed cavity with inhomogeneity in the z-direction. 
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